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Abstract 



We give a semiclassical interpretation of the mass asymmetry in the fission of heavy nuclei. 



Using only a few classical periodic orbits and a cavity model for the nuclear mean field, we 
reproduce the onset of left-right asymmetric shapes at the fission isomer minimum and the 
correct topology of the deformation energy surface in the region of the outer fission barrier. We 
point at the correspondence of the single-particle quantum states responsible for the asymmetry 
\ with the leading classical orbits, both lying in similar equatorial planes perpendicular to the 
symmetry axis of the system. 
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One characteristic feature of the fission of many actinide nuclei is the asymmetric mass 
distribution of the fission fragments. The liquid drop model 0, although able to describe 
many aspects of the fission process qualitatively, cannot explain this mass asymmetry in heavy 
nuclei where the fissility parameter x is close to unity 0: the balance between the attractive 
surface tension and the repulsive Coulomb force favors left-right symmetric shapes and thus also 
the symmetric fission. An explanation for the observed asymmetry of the fragment distributions 
became possible with Strutinsky's shell-correction method || which includes the quantal shell 
effects stemming from the discrete spectra of the nucleons in their mean fields. One writes the 
total binding energy of a nucleus consisting of N neutrons and Z protons as 

E tot (N,Z;def) = E LDM (N, Z; def) 

+5E n (N;def) + 5E p (Z;def). (1) 

Here E LDM is the liquid drop model (LDM) energy; 8E n and 5E p are the shell-correction 
energies of the neutrons and protons, respectively, which are obtained in terms of the single- 
particle energies of realistic shell-model potentials. All ingredients depend on the shape of the 
nucleus, which is described by some suitable deformation parameters, summarized in ([!]) by 
'<ie/'. The shell-correction approach was very successful in reproducing experimental nuclear 
binding energies and fission barriers @, ||, |6| at times where purely microscopical selfconsistent 
calculations of Hartree-Fock type were not yet available |7j . 

In Fig. |l] we show a schematic fission barrier of a typical actinide nucleus, taken along the 
adiabatic path through the multi-dimensional deformation space. The heavy dashed line is the 
LDM deformation energy which leads to a spherical ground state and to equal fragment sizes. 
The solid line is the total energy according to Eq. ([I]); it contains the typical shell oscillations 
coming from the shell corrections SE n and SE p . These lead to a deformed ground state minimum 
and to a higher-lying local minimum, the so-called fission isomer. (For an extended review on 
the physics of the 'double-humped fission barrier', see Ref. ||.) The shapes assumed hereby 
have axial symmetry and left-right symmetry. When the latter is relaxed, the energy is found 



||, 10] to be lowered along the way over the outer barrier, starting at the fission isomer. The 
gain in energy persists all the way down towards the scission point, where the nucleus breaks 
into two fragments of unequal size. Nonaxial deformations do not change this feature; they 



only lead to a slight reduction of the inner barrier (see Ref. JTTJ for a short review of fission 
barrier calculations and the role of various deformations). It is important to note that the 
onset of the mass asymmetry takes place already at an early stage of the fission process, long 
before the nucleus breaks up. It is a pure quantum effect which only comes about if the shell 
corrections are included into the total energy. The microscopic origin of this instability has 



been investigated by Gustafsson, Moller, and Nilsson [0. They found that only two specific 
types of single-particle states with large angular momenta along the symmetry axis are strongly 
sensitive to the left-right asymmetry: one of them has the maxima of its wave functions along 
the central waist line of the nucleus (see the upper right in Fig. [I]), whereas the other has 
maxima along the circumferences of two equatorial planes at some distance of the center (with 
opposite phases on either side, see the middle right in Fig. [TJ). The coupling of these states 
through the left-right asymmetric components of the mean field leads to a decrease of one set 
of eigenenergies which lie below the Fermi energy, and thus to a reduction of the total binding 
energy when the asymmetry is switched on (see the lower right in Fig. [I]). 

In the following we give a very simple and transparent semiclassical interpretation of this 
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Figure 1: Left side: schematic fission barrier of a typical actinide nucleus (after Ref. JTTJ). Right 
side: schematic probability amplitudes of the leading single-particle wavefunctions responsible 
for the onset of mass asymmetry (after Ref. fll2"|). 



quantum effect. We employ the periodic orbit theory (POT) developed independently by 
Gutzwiller ]TJ| for systems with isolated classical orbits and by Balian and Bloch JTJJ for 
arbitrarily deformed cavities with reflecting walls. In the POT, one obtains the oscillating part 
of the level density of a quantum system from the so-called 'trace formula': 



8g (E)~J2Ao(E) cos 



po 



Spo(E) 



71 



a 



po , 



(2) 



The sum is taken over all periodic orbits, labeled 'po', of the corresponding classical system. 
S po are the classical actions along the periodic orbits and a po are phases related to the number 
of conjugate points along the orbits. The amplitude A po of each orbit depends on its period, its 
stability and its degeneracy. Together with the smooth part g (E) which can be obtained in the 
(extended) Thomas-Fermi model, it approximates the exact quantum-mechanical level density: 
g (E) + 5g (E) = g (E) = J2i 8{E — Ei), where E^ are the eigenenergies of the system and the 
sum runs over all quantum states i. (See Ref. |D| for an introduction to the POT and detailed 
explanations of all the above ingredients). Gutzwiller's trace formula has initiated promising 
steps towards the semiclassical quantization of chaotic systems |l(J. A different use of the POT 
15, O] is to obtain a coarse-grained level density by keeping only the shortest orbits with the 



largest amplitudes in the trace formula (0). This allows one to relate the gross-shell structure 
of interacting fermion systems in the mean-field approximation to a few classical orbits. Using 
Eq. 



the semiclassical expression for the energy shell correction 5E becomes JT7] 
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Here E F is the Fermi energy and T po = dS po / dE\ Ep the period of the orbit labeled po. Pioneering 
work has been done in this direction by Strutinsky et al. [0, who generalized the Gutzwiller 
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theory to systems with continuous symmetries and used it to give a semiclassical explanation 
of the systematics of nuclear ground-state deformations. They showed that the ground-state 
minima of the nuclear binding energy as functions of particle number and deformation follow 
the loci of constant action of the shortest periodic orbits. Another beautiful example is the 
beating pattern of the coarse-grained level density in a spherical cavity, which was related by 
Balian and Bloch [14] to the interference of the triangular and square periodic orbits, and later 
predicted |Tj| and observed |B| in metal clusters in the form of the so-called supershells. 

For our semiclassical investigation, we replace the nuclear mean field by a cavity with 
reflecting walls and consider only one kind of particles. This should yield the main physical 
effects, since the neutron contribution 5E n to the total energy contains the largest part of the 
shell oscillations, as seen, e.g., in Fig. VIII-4 of Ref. ||]. We employ the parameterization 
(c, h, a) of this reference, used there to define the shape of the liquid drop. Here it defines the 
boundary of the cavity in cylindrical coordinates (p, z, 0), with the symmetry axis along z, by 
the shape function p = p{z\ c, h,a). c is the length of the semiaxis along the z direction in units 
of the radius Rq of the spherical cavity (c=l, h=a=0). h is a 'neck parameter' that regulates 
the formation of a neck leading to the scission of the nucleus into two fragments. a^O yields 
left-right asymmetric shapes. The volume of the cavity is kept constant for all deformations. 
(See Ref. |3J] for the details of this parameterization, and especially Figs. VII- 1 and VIII-5 for 
the most important shapes occurring in fission.) The parameters (c, h) are chosen such that the 
one-dimensional curve h=a=0 along c follows the adiabatic fission barrier of the LDM (shown 
schematically in Fig. |T]). Even including the shell effects, h=0 gives a reasonable picture of the 
double-humped fission barrier. 

We now have to determine the shortest periodic orbits of this system to calculate the 
gross-shell structure in 5E. At large deformations (here c ~ 1.4), these are the orbits lying 



in equatorial planes perpendicular to the symmetry axis |20], pi] . The positions Zi of these 
planes along the z axis are given by the condition that the shape function be stationary: 
dp(z; c, h, a)/dz\ Zi = 0. The periodic orbits have the form of regular polygons and are char- 
acterized by (p, t), where p is the number of reflections at the boundary and t the number of 
windings around the symmetry axis. Both numbers are restricted by p > 2 and t < p/2. The 
semiclassical contribution of such orbits to the trace formula (0) has been derived by Balian 
and Bloch |TJ| ; we refer to their paper for the explicit form of the amplitudes Apt and phases 
a pt . The lengths of the orbits are L p } = 2pRi sm(irt/p), where Ri = p(zf, c, h, a), so that their 
actions equal Sp^(E F ) = hk F L p ^ in terms of the Fermi wave number k F = y/2mEp/h. 

The range of validity of Eqs. (0) and (B|) is, however, limited. They are correct only as long 
as the orbits are well separated from neighboring periodic orbits, in particular as long as the 
orbits are not close to a bifurcation. At a bifurcation the amplitudes A pt diverge and the trace 
formula has to be modified. Generally, bifurcations exist in different forms, but for the shapes 
studied here we need consider only one type of bifurcation. It occurs when the positions z { of 
several equatorial planes coincide. In the (c, h, a) parameterization, there are at most three 
such planes. One plane always exists; the other two arise at the points (co, ho, ao) where the 
neck formation starts. In the symmetric case (a=0), one plane is always located at zq=0 and, 
beyond the bifurcation point, the other two are located symmetrically at rkz\ (with zi>0) and 
contain identical periodic orbits. 

Near a bifurcation point, the neighboring orbits give a joint contribution to the level density. 
This contribution can be derived from an oscillatory integral which contains the contributions 
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of all orbits of type (p,t) from the different planes. The integral is of the form [14 



dz fp t (z) exp{ik F L pt (z)} 



(4) 



where f p t(z) are slowly varying analytic functions of z. Since the plane positions Zi of the peri- 
odic orbits are determined by the stationary points of the length function L pt (z) = 2p sin(7rt/p) • 
p(z;c,h,a), a stationary phase evaluation of leads back to separate contributions to Eq. 
(3) for each plane, with the amplitudes and phases given in [14|. In order to obtain an ap- 



proximation to (H) that is valid at the bifurcation as well as far from it, we employ a uniform 
approximation that is appropriate for the case of three nearly coincident stationary points in a 
one-dimensional oscillatory integral [22] . It is expressed in terms of Pearcey's integral and its 
derivatives: 



SE pt = 3fte {[u A P(ui,u 2 ) + u 5 P x (u 1 ,u 2 ) 
+u 6 P y ( Ul ,u 2 )}e iU3 } , 

where Pearcey's integral is defined by 



P(x,y) 



dz exp[i(z 4 + xz 2 + yz)], 



(5) 



(6) 



and P x (x,y), P y (x,y) denotes its derivatives with respect to its first and second argument, 
respectively. The real constants u\ . . . Uq are determined by the semiclassical amplitudes, ac- 
tions and phases of the periodic orbits. If the orbits are well separated, Eq. (H) reduces to 
contributions to the standard trace formula fl3|). 

When the asymmetry parameter a is sufficiently large, the contributions from one plane of 
orbits can be evaluated in the stationary phase approximation. Both in this limit and in the 
symmetric case (a=0), the result can be expressed in terms of cylindrical Bessel functions J^x), 
and the formulae are analogous to those for generic saddle-node and pitchfork bifurcations 
We give here the result for 5E pt in the symmetric case: 
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Here L 



[1$ + 4° } ]/2 and AL pt = [L W 



vt 



Lp t ]/2 in terms of the lengths 1$ , Lp L t ' of the 
orbits pt situated at z = z and z = ±2i, respectively, v equals —1 before the bifurcation (i.e., 
for only one orbit plane) and +1 after the bifurcation (for three orbit planes). 

In the right-hand panels of Fig. [|we show contour plots of the semiclassical shell-correction 
energy SE in the (c, a) plane for two values of the neck parameter. The energy unit is Eq = 
h 2 /2mRQ, where Rq is the radius of the spherical box. The Fermi wave number k F was chosen 
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Figure 2: Contour plots of deformation energy versus elongation c and asymmetry a for two val- 
ues of the neck parameter: h=0 (above) and h=— 0.075 (below). Leftside: quantum-mechanical 
neutron shell-correction energy of 240 Pu from Ref. H. Right side: semiclassical shell-correction 
energy SE. White dashed lines indicate the loci of constant classical action of the central 
equatorial periodic orbits; black dashed lines the loci of the bifurcation points. 

such that SE has a minimum at the deformation c =1.42, h=a=0 of the fission isomer. On the 
left of Fig. ^| we have reproduced the neutron shell-correction energy 5E n of the nucleus 240 Pu, 
obtained in Ref. || with a realistic Woods-Saxon type shell-model potential. We see that the 
semiclassical result correctly reproduces the topology of the deformation energy in the (c, a) 
plane for both values of h, in particular the onset of the mass asymmetry at the fission isomer. 
It should be noted that we have only included orbits with winding number one (t=l) and with 
up to Pmax = 10 reflections. The results for SE are the same within a few percent when only 
orbits with p = 2 and 3, i.e., only the diameter and triangle orbits, are included. The loci of 
the bifurcation points (c , a ) are indicated in Fig. |2] by the black heavy dashed lines (that for 
h=— 0.075 is hardly visible in the upper right corner of the plot). This shows that the essential 
feature, namely the energy gain due to the asymmetric deformations for c > 1.42, is brought 
about by only two classical orbits: the diameter and the triangle in the central equatorial plane. 
The white dashed lines give the loci of constant actions of the periodic orbits at zq, fixing their 
value at a=0. (Note that the actions of all orbits in a given equatorial plane have the same 
deformation dependence.) We see that the valley that leads from the isomer minimum over 
the outer fission barrier in the energetically most favorable way is following exactly the path 
of constant action of the leading classical orbits; the path is practically identical with that 
obtained in the quantum-mechanical shell-correction calculations. A more detailed comparison 
with the latter, together with the calculational details of our studies, will be presented in a 
forthcoming publication. 
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Figure 3: Right side: semiclassical shell-correction energy SE as in Fig. for h = 0, but in a 
three-dimensional perspective view. Arrows 'symmetric' and 'asymmetric' show two alternative 
fission paths. Left side: shapes found along the asymmetric fission path; the equatorial planes 
of periodic orbits included in the trace formula are shown by the vertical lines (solid lines for 
stable and dashed lines for unstable orbits). 



In Fig. [3], we show the same results as in the upper right part of Fig. ^ but in a perspective 
view of a three-dimensional energy surface. On the left, the shapes of the cavity are given 
corresponding to the points A at the isomer minimum and the points B and C along the 
asymmetric fission barrier (see also the corresponding points in Fig. 0). Note that C lies 
beyond the bifurcation point and thus contains three planes of periodic orbits. 

In summary, we have shown how a specific quantum effect, causing a drastic rearrangement 
of the shape of a complex many body system, can be described semiclassically by the constancy 
of the actions of a few periodic orbits. We point out the close correspondence of the equatorial 
planes of these orbits (see the left of Fig. |3|) with the locations of the wave function maxima (see 
the right of Fig. [I]) of the relevant quantum states. We also note that the classical dynamics 
of this system is rather chaotic, particularly in the asymmetric case and for small values of the 
conserved angular momentum L z [p4] . The regular regions in phase space, connected to the 



stable periodic orbits with axial degeneracy, are still strong enough to cause this important 
shell effect. 
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